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Abstract 

To study the shape formation process of carbon nanotubes, a string equation 

describing the possible existing shapes of the axis-curve of multishell carbon 

tubes (MCTs) is obtained in the continuum limit by minimizing the shape 

energy, that is the difference between the MCT energy and the energy of 

the carbonaceous mesophase (CM). It is shown that there exists a threshold 

relation of the outmost and inmost radii, that gives a parameter regime in 

which a straight MCT will be bent or twisted. Among the deformed shapes, 

the regular coiled MCTs are shown being one of the solutions of the string 

equation. In particular, the optimal ratio of pitch p and radius ro for such 

a coil is found to be equal to 2tt, which is in good agreement with recent 

observation of coil formation in MCTs by Zhang et al. 
61.48.+C, 68.65.+g, 68.70.+W 
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Since the discovery of straight and multishell carbon nanotubes (MCTs) in arc discharges 
|jlj], many unique and novel properties have been predicted for the tubes. Among them, an 
especially intriguing one is their structural stability, the mechanical properties of MCTs are 
expected to be significantly stiffer than any presently known materials |2|. However, in the 
recent synthesis by the catalytic decomposition of gas such as acetylene, a significant fraction 
of the produced MCTs exhibits various curved shapes |§, of which the most striking shapes 
are helices, i. e., the regular coils. It has been pointed out that the tubes can be twisted 
and deformed by an abrupt release of energy and a singularity in the stress-strain curve [§J . 
In Ref. ||, the regular coil formation was explained by a periodic distribution of pentagon- 
heptagon-pair dislocations (PHPDs), but why the PHPDs are distributed periodicly rather 
randomly still remains as an open question. Dunlap pointed out that the regularity could 
be caused by defect-defect interactions and constraints on the optimal heptagon-pentagon 
nanotube bend ||. In Ref. 0, the bent carbon tubes were simulated by classical molecular 
dynamics based upon the three-body Tersoff-Brenner interatomic potential, but it is rather 
difficult to extract the essential physics from such a numerical approach. 

Therefore, the general questions can be posed as follows, what is the mechanism of the 
curved deformation for the MCTs, and can we derive the deformed shape. 

In this Letter, we analytically obtain the general equilibrium-shape equation of the axis- 
curve of the MCT in the continuum limit by taking account of competition among the 
curvature elasticity, the adhesion of the interlayer van der Waals bonding and the tension of 
the outer and inner surfaces of a MCT. The sum of these three energies can be understood as 
the shape formation energy (see bellow for details). We find that the variation of the shape 
formation energy yields an equation of rigid string, which has been studied extensively 
in differential geometry 0. Straight line and regular coil are two exact solutions of the 
equation. We show that under certain geometric conditions the shape formation energy of a 
straight MCT could become negative, in other words, the straight MCT becomes unstable 
in a quench-like formation process, and as a result, bent or twisted MCTs will be formed 
spontaneously to keep the equilibrium condition (i. e. the zero shape formation energy). 
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Taking into consideration of the equilibrium condition in the quench-like cooling processes, 
the above argument provides an insight for the mechanism of the curved deformation for 
the MCTs. The optimal ratio of the pitches and radii of the regular coils formed in such 
processes is shown impressively to be equal to 2n. Our result is in good agreement with 
recent observation of the coil formation in MCTs by Zhang et al. [|[]. 

Generally, when the hydrocarbons are thermally decomposited |§ , the carbon molecules 
are condensed mostly to form an isotropic smectic-like crystal, i. e. the carbonaceous 
mesophase (CM), while the remained space is filled by plate-like molecules f|. Thus the 
MCT formation is quite similar to the tube formation of a smectic-A phase grown from 
isotropic phase in liquid crystal ||. As shown in Ref. ||, the shape formation energy is 
the additional energy of an MCT with respect to CM, which is a sum of the following 
three terms, (i) the net difference of the volume free energy between MCT and CM, i.e. 
Fy = —goV where V is the volume of the MCT and —go is the adhesion energy density of 
the interlayer van der Waals bonding, (ii) the surface energy Fa = 7(v4 + v4j) where 7 is the 
surface tension, A a and Ai are the areas of the outmost and inmost surfaces, respectively, 
and (iii) the curvature elastic energy of the layers. 

We consider first the third term of the shape formation energy, the curvature elastic 
energy. MCTs can be treated as a set of curved graphite layers 0. For a single layer, the 
curvature elastic energy is an increment part of the in-layer covalent energy due to the layer 



curvature. Following Lenosky et al. jTOj? the curvature elastic energy of a single layer curved 
graphite carbon has the form as, 

E b S) = e l ZXZ! U v) 2 + C 2 H I 1 - n i • n i) + e 3 ( n i • U ij)( n i • U Ji), W 

i <j> <i,j> <i,j> 

where Ujj is a unit vector pointing from carbon atom i to its neighbour j and rij is a unit 
vector normal to the fullerene surface at atom i. The summation J2<j> is taken over the 
three nearest neighbour j atoms to atom i, and the sums of the last terms are taken over 
only the nearest neighbour atoms. The superscript (s) emphasizes the energy Ef, is for a 
single layer. Our first task is to reduce Eq. (1) into a continuum form. 



A curved single-shell tube of radius p without inclusion of its two end-caps can be de- 
scribed by 

Y(s,<f>) = r(s) + p(N(s)cos0 + b(s)sin0), (2) 

where < <p < 2tt, and < s < I is the arc- length parameter along the curved tube- axis 
whereas curve r(s) is the vector representation for the curve of the tube axis. I is the total 
length of the tube. N(s) and b(s) are the unit normal and unit binormal vectors of r(s) 
respectively Making use of the well-known Frenet formulas [fQ |, 



t a = k(s)N, N s = -Jfc(s)t - r(s)b, b s = r(s)N, (3) 

where t = dr/ds, t s = dt/ds, N s = dN/ds, h s = dh/ds, k(s) and r(s) are the curvature 
and torsion of r(s) respectively. In ref. [0, we have derived the area element cL4 = p(l — 
pk cos <j))d<pds, the mean curvature = (2pkcos<p — l)/2p(l — pk cos 4>) and the Gaussian 
curvature = — fccos0/p(l — pk cos <p) for the tube surface of Y. It is then easy to prove 
§ KdA = 0. According to the Gauss-Bonnet theorem [IT], this means the topology of the 



curved tube is the same as that of a straight tube. From the Euler's theorem in topology, 
the surface can be perfectly embedded by a carbon network of six-member rings as in plane 
graphite layer. With the help of Frenet formula (3), we transform the vector functions in Eq. 
(1) into continuum limit by expanding them up to the order of 0(a 2 k 2 ), where a = 1.42A is 
the bond length in the graphite layer, 

Ui (M) = Uij = (1 - ^k 2 {M))t{M) + £k(M) + ^fc s (M))N(M) - ^k(M)r(M)b(M), 
o 2 o o 

(4) 

where k s = dk/ds, M — 1, 2, 3 denote three families of sp 2 -bonded curves with one curve of 
each family acrossing from atom i to its three neighbour atoms j on the surface where the 
carbon atoms embedded, and have one to one correspondence to atoms j. Here we would like 
to emphasize that all the expressions K(M), t(M) ,N(M), r(M), and b(M) are functions 
of the arc- length s, where s = ia. Specifying the tube surface described by Eq. (2), in 



which the sp 2 -bonded curves can be considered approximately as geodesies for the curved 
surface Y, we have additionally N(M) = n i; k(M) = c x cos 2 6{M) + c 2 sin 2 6 l (M), r(M) = 
(ci — c 2 ) sin 6(M) cos 6(M) where c\ and c 2 are the two principal curvatures of the surface 
at atom i location, i. e., H = (c± + c 2 )/2 and K = cic 2 , and 9(M) are the angles between 
ci direction and t(M). Considering ELi sin2 #( M ) = 3 / 2 and Em=i sin4 #(M) = 9/8, and 
substituting Eq. (4) into Eq. (1), we obtain an important formula for the curvature elastic 
energy of the tube pf, 



E ( b s) = j\\um 2 + kK]dA , (5) 

where the bending elastic constant, 

k c = (1/32) (18 £l + 24e 2 + 9e 3 )(a 2 / ( x) , (6) 

with a = v / 3« 2 /4 = 2.62A 2 being the occupied area per atom, and the saddle-splay modulus 
is. 

k = -(8e 2 + 2e 3 )k c /(Qe 1 + 8e 2 + 3e 3 ) . (7) 

Formula (5) is actually a general expression of the elastic energy which is valid also to fluid 
membranes f[3| and solid shells fL4|. If we substitute (ei,e 2 ,es) in Eqs. (6-7) by the values 



of (0.96, 1.29, 0.05) eV respectively, which were calculated by Lenosky et al. using a local 
density approximation |]10"| |, we find k c = 1.17 eV and k c /k = —1.56. The obtained value of 



k c is in reasonable agreement with the value of 1.02 eV calculated by Tersoff |L5| using an 
atomistic method for straight tubes, and is excellently close to the value of 1..2 eV extracted 
from the measured phonon spectrum of graphite [16|]. The calculated ratio of k c /k is also 



close to the result of k c /k = —105.4/88 = —1.2 measured by Blakeslee et al. [17]. Therefore, 
we have sufficient confidence on Eq. (5). Moreover, since it has been averaged over three 



nearest neighbours for each site in expression Eq. (1) of Lenosky et al. |TIJ, we can have only 
two invariants as H 2 dA and KdA in Eq. (5). Therefore, within the same approximation, 
the free energy corresponding to Eq. (5) should be again a linear combination of these 



two invariants only. Consequently, would have the same formal expression as Eq. (5) 
with coefficient k c and k being temperature dependent. 

To extend the above result to a MCT, one has to integrate Eq. (5) from its inmost radius 



Pi to the outmost radius p Q . We may apply a similar treatment as that of [18] which is in 
fact devoted to the curved smectic crystal multilayers, and has a layer structure similar to 
graphites'. Replacing fc n with k c /d in Eq. (3) of Ref. and neglecting the constant term 
associated with j> KdA, we have 

F b = £if> = (k c /2d) / 2Vm=KM \~_ D D H H + _ D D V J^ )dA, (8) 

where d = 3AA is the space between two neighbour graphite layers and D = p a — pi is the 
thickness of the MCT. Here, the surface integral is carried out over the inner surfaces. Using 
the above expressions for H, K, and dA and integrating from = to 2tt, we obtained the 
curvature elastic energy for the MCT as 



{-Kkc/d) 



ds. (9) 



Pi l + ^l-k 2 pl 

We now turn to consider the other two terms of the shape formation energy, Fy and 
Fa, both of which are weak binding energy, i. e. the adhesion energy between layers of an 
MCT. Despite the fact that many of the structural properties of plane graphites are well 
understood, the calculation of interlayer adhesion energy for curved graphites is still an open 
question. The observation in Ref. pj reveals that the interlayer distance d in MCTs remains 
to be the same as that in plane graphite, but the in-layer lattice structures for each single- 
shell tube in one MCT may have different helicity. In other words, the interlayer lattices are 
not in perfect registry (referred to "incommensurate" or "mismatched" lattices). Therefore, 
the attractive forces between layers in an MCT cannot be accounted for by conventional 
forces in plane graphites |19| . However, the adhesion energy for mismatched lattices is often 
smaller than that for commensurate surfaces. We use a mismatched parameter rj to take 
account for the mismatched effect between the interlayer lattices, < r\ < 1, and rj — 1 
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corresponds to the commensurate case. Since we are so far not aware of more detailed 
knowledges, as the lowest approximation, we take the following simple energy form, 

F V + F A = -g 7i(p 2 - pi) J ds + 27T 7 (po + Pi) J ds, (10) 

where —go = t]AE c /d, AE C = —330erg/cm 2 = —2.04 eV/nm 2 is the interlayer cohesion 
energy of planar graphite obtained theoretically by Girifalco and Lad and 7 is the 
tension of the outmost and inmost surfaces of the MCT, which is equal to half of the energy 
needed to separate two unit-area surfaces, i. e., 7 = — (1/2) A^. 



Usually we have the MCTs with p 2 Q k 2 « 1, then the expression ln[(l + y 1 — k 2 pf) /(l + 
1 — k 2 p 2 )] can be approximated to (1/4) (p 2 — p 2 )k 2 , and the shape formation energy of 
the MCT, Eqs. (9-10), can be subsequently simplified to, 

F = F v + F A + F b = mJds + aJ k 2 ds , (11) 

where m = ir(k c /d) \n(p Q / p { ) + 27T7(p + p{) - 7rg (p 2 o - pf) and a = (1/ '4)n(k c / 'd)(p 2 - pf). 
Eq. (11) is nothing but a string action 2l|. The variational equation 5F = yields the 
equilibrium-shape equations of the string 0, 



m . 



2k ss + k 3 - 2kr 2 k = 0, (12) 



k 2 r = const., (13) 

where k ss = d 2 k(s)/ds 2 . 

Following what has been discussed for the derivation of curvature elastic energy, in the 
derivation for these three terms of the shape formation energy, we have expressed all the 
relevant quantities in terms of geometric language. Therefore, due to the reason of geometric 
symmetry on the low dimensional manifolds, the corresponding free energy should have 
the same expression as those of the derived shape formation energy with the coefficients 
becoming temperature dependent. We will keep such understand in the above and hereafter 
discussions. 



It is obvious that a straight line is always a solution of the string equations (12) and (13), 
since its k and r are zero. The corresponding shape formation energy of a straight MCT 
is F = ml. The shape formation energy is regarded as a free energy and the equilibrium 
threshold condition of F = yields the criteria for the growth of a straight MCT as 

m = n(k c /d) \a(p / p{) + 2^(p + p { ) - irg {pl - p- ) = . (14) 

This equation describes the geometric relation between p Q and pi in terms of the physical 
parameters k c , 7, and go for a straight MCT. Both 7 and go are also dependent on the 
formation temperatures and catalyst. So the detailed data measured from the produced 
MCTs can reveal the properties of 7 and go with the help of Eq. (14). The formation 
procedure for MCTs, either a quick growth in which the temperature can be regarded as 
constant or a sudden cooling, is actually a sort of quench-like process. As long as the 
shape formation energy for the straight MCTs, i. e., deviates downwards from the threshold 
condition in the formation procedure, becomes negative, the resultant remnant part of energy 
will prevent the straight MCT from keeping stable. Then a shape deformation will be 
induced and it would lead to another solution of the string equation with its shape formation 
again being equal to zero. Therefore, by considering the threshold condition F = 0, any 
outward growth by increasing po will make the straight MCT undergo a shape deformation as 
long as go being kept constant. Furthermore, go may increase with temperature decreasing, 
following again Eq. (14), a straight MCT grown may also be coiled under the cooling process. 
These features give a natural explanation for the deformation of MCTs. 

Now we would prove that the observed MCT regular coils shown in Ref. are just 
the allowed solutions of equations (12) and (13). Mathematically, the regular coils can be 
described by vectors, 

r(s) = (r cosus, r sinus, hus), (15) 

where the coiled pitch p = 2nh, ro is the coil radius and R = oj~ 1 = \J r\ + h 2 . One can 
easily show from the Frenet formulas, Eq. (3), that k = w 2 ro, r = —u 2 h, and the regular 
coil curves are the solutions of (12) and (13) if their r and h satisfy the following equation: 
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r 2_ 2 /i 2 --(r 2 + / i 2 ) 2 = 0. (16) 



a 



Introducing h/R = sin8, and r /R = cos#, and taking into account of Eq. (16), we have 
R 2 = (ct/m)(cos 2 #-2sin 2 #), and k 2 = r^/R 4 = (m/a) [1/(1 - 2 tan 2 6)}. Therefore, the coil 
formation energy can be derived from (11) as 

F = m V + T=Lri- (l7) 

Since now we have the coil situation which is quite different from the straight MCT case, 
even for the negative value of m in Eq. (17), we may treat the threshold condition F = 
as tan 2 6 = h 2 jr\ = 1 or 

^ = 2tt . (18) 
r 

We compare the optimal ratio given by Eq. (18) with the experimental results reported 
in Ref. and find a good agreement. As shown in Fig. 1 of Ref. J|, there is a fraction 
(about 10%) of MCTs being regularly coiled with a variety of radii tq and helix pitches 
p. By a direct evaluation from the figure, we do find p/r = 2tt is well hold. A rough 
estimation from the coil shown in the inset of this fig 600 nm and r$ « 100 

nm, i. e. p/ro ~ 6 ~ 2ir. Another coil, shown in Fig. 2 of the same reference, has its 
p ~ 700 nm and r ~ 100 nm which leads to p/r w 7, again close to the prediction of 
Eq. (18). Moreover, the results corresponding to the coil are shown in Fig. 3 of Ref. 
tan 9 = f — 0o ~ 2.2/2 w 1.1, are also in good agreement with the present prediction of 
tan# = 1, where 0o is defined within the context of Ref. ||. 

It is also interesting to study the value of m for the MCT coil by utilizing the data 
provided in Ref. 0, where the experimentally observed values of 2pi = 3 ~ 7 nm and 
2p a = 15 ~ 20 nm. Making use of calculated values, k c = 1.17 eV and g = rj x 2.04 x 
10~ 2 (eV/A 2 )/d = 6.0377 eV/nm 3 , we can estimate ng (pl — p 2 ) = (I.O77 ~ IJ77) xl0 3 eV/nm, 
and 7r(k c /d)ln(p / pi) + 2^{p Q + p { ) = 68 ~ 101 eV/nm. Considering expression (14), we 
find that under reasonable approximation, m in the practically formed MCT coils always 



9 



take negative value. Such a fact is notably consistent with the above proposed explanation 
for the shape deformation mechanism of M CTs. 

In summary, by deriving a string action type expression for the formation energy of the 
MCTs as well as its equilibrium-shape equation, we have shown that there is a threshold 
condition for the formation of straight MCTs, below that the straight MCTs become unstable 
and will undergo a shape deformation. In particular, we derive further an optimal formation 
condition p/r = 2n for the regular coil solution which is in good agreement with the recent 
experiment observations. 

This work is partly supported by the National Natural Science Foundation of China. 
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